Vanishing conductivity of quantum solitons in polyacetylene 
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Abstract 

Quantum solitons or polarons are supposed to play a crucial role in the electric conductivity of 
polyacetylene, in the intermediate doping regime. We present an exact fully quantized calculation 
of the quantum soliton conductivity in polyacetylene and show that it vanishes exactly. This is 
obtained by applying a general method of soliton quantization, based on order-disorder duality, to 
a Z(2)-symmetric complex extension of the TLM dimerization effective field theory. We show that, 
in this theory, polyacetylene solitons are sine-Gordon solitons in the phase of the complex field. 

PACS numbers: ll.10.Kk, ll.10.Lm, 61.82.Pv 



*Electronic address: leo.mondaini@pq.cnpq.br 
^Electronic address: marino@if.ufrj.br 
^Electronic address: alex@lana.ccne.ufsm.br 



1 



I. INTRODUCTION 



The discovery of a tremendous increase in the electrical conductivity of trans- 



polyacetylene, when doped either with halogens or alkalis, [1[ was a breakthrough of far 
reaching consequences in physics and chemistry. The fact that the trans-isomer occurs in 
two degenerate species opens the possibility of occurrence of soliton defects interconnecting 
them. This fact unfolded an enormous range of possibilities interconnecting many areas, 
including mathematics, theoretical and experimental physics. The subject has remained on 
the focus of interest until recently [2J. It actually happens that such topological excitations 
are produced in the process of doping 0, 4\. Pure polyacetylene has one active 7r-electron 
per site and is a Peierls insulator, due to the electron-lattice interaction. It has been found 
that in the presence of a soliton, treated at the classical level, electron states are created in 
the middle of the gap, hence it is energetically favorable for the extra doped electrons to 
create solitons and occupy the midgap states rather than going into the conduction band. 

From the very beginning the existence of three different doping regimes became clear. 
Firstly, for low doping concentrations (small compared to 1%) the above picture of classical 
solitons works very well. The solitons are pinned by the dopant atoms, which create them 
with the corresponding midgap states. The conductivity is thermally activated, correspond- 
ing to a transition from the midgap to the conduction band states and can be understood 
quite similarly to the conductivity in semiconductors. Secondly, for high concentration of 
dopants (5% to 10%) the conduction regime is clearly metallic, with an unfilled conduction 
band, and can be thereby understood. Thirdly, there is an intermediate regime of dopant 
concentration, of the order of 1%, in which none of the previous models work. In this regime, 
the solitons become dynamic carriers of charge and a full quantum treatment of these exci- 
tations becomes unavoidable since the soliton mass is of the same order of the electron mass 

I 

The purpose of this work is to apply a general method of quantization of soliton excita- 
tions [6], in order to describe the conductivity of polyacetylene in the intermediate regime. In 
order to do that, however, the following obstacle must be removed. From the mathematical 
point of view, the soliton is a topologically nontrivial configuration of the dimerization - or 
phonon - field, which describes the lattice degrees of freedom. This, of course, is a real field, 
whose effective potential has Z(2)-symmetry and two degenerate minima that correspond to 
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the two species of trans-polyacetylene. The above method of soliton quantization, however, 
only applies to complex fields, in the case of a multiplicative symmetry such as Z(2) 6|. 

Therefore, in order to describe the quantum solitons of the system and specifically their 
role in the electric conduction in the intermediate doping regime, we propose a Z(2)- 
symmetric complex extension of the effective potential for the dimerization field. This has 
the same topological properties of the former and, consequently should not alter substan- 
tially the soliton physics. 

In Section 2, we describe the method of soliton quantization in a theory of a complex 
scalar field with Z(N) symmetry and show that the soliton excitations of it are sine-Gordon 
(SG) solitons in the phase of the complex scalar field. In Section 3, we propose the N = 2 
version of this theory as the complex extension of the effective theory for the dimerization 
field in polyacetylene. In Section 4, we derive an exact series expression for the quantum 
soliton current- current correlation function and, out of it, obtain the soliton conductivity. 
We show that this exactly vanishes. In conclusion, we have an exact demonstration at a full 
quantum level, that dynamic solitons are actually not the carriers of charge in polyacetylene. 
Rather, polarons, which are basically soliton-antisoliton bound states should be responsible 
for the transport of charge in the intermediate doping regime of polyacetylene. 



II. QUANTUM PHASE SOLITONS IN THEORIES WITH A Z(N) SYMMETRY 

We start by considering the following theory describing a complex scalar field in (1+1)- 
dimensions, 

£ = d^Ptf, + 7 {<P* N + <P N )~V (0*0) M , (1) 
where N and M are integers and 7 and rj are real parameters. This is invariant under the 

■ 2-7T 

Z(N) transformation: (fi(x,t) — ► e % ~^ (f>(x, t). The choice 7 > 0, implies the spontaneous 
breakdown of the Z(N) symmetry. In this case, the theory will have degenerate vacua and 

nnn 

soliton excitations. A full quantum theory of these solitons was developed in [6|, |7|, lSj. This 
includes an explicit expression for the soliton creation operator, namely 



= exp |~ £Ve M V(0i<K0} , (2) 



and a general expression for its local Euclidean correlation functions 



,13, 



( i i{x) t i\y))=N J P0*P0expj- j tPzKD^YiD^ + V (</>*,. 



(3) 



where 

D^ = d^- iaA„, A„(z, C) = [" d^ v 8\z - £). (4) 

Jx,C 

In the above expression, V is the potential of an arbitrary Lagrangian and the integral is 
taken along an arbitrary curve C, connecting x and y. It can be shown, however, that (J3]) 
is independent of the chosen curve. 

We are going to show in what follows that these quantum solitons may be identified with 
SG quantum solitons in the phase of the field 0. 

Using the polar representation for 0, namely, 4>(x,t) = p(x,t)e ( x,t \ where p and 9 are 
real fields, we can rewrite the above Lagrangian as 

£ = d^p + p 2 d ll 9d»9 + 2 1P N cos N9 - 7] P 2M . (5) 

In what follows, we will be interested in the topological properties of the theory. As we 
shall argue, these are not affected by p fluctuations, hence, from now on we will make the 
constant p approximation, 

p(x,t) = p , p constant. (6) 

Using this in ([5]) we get 

C = pld^O + 2 7 p Ar cos N9 - V pl M , (7) 

which is a SG Lagrangian in 9. 

We conclude that, in the constant-p approximation, the theories given by ([1]) will present 
SG solitons in the phase of the complex scalar field 0. The corresponding topological current 
will be 

J» = e^OJ, (8) 
which is associated to the topological charge operator 

/oo /»oo 
dx'J°= dx'd x ,e(x',t) = 9(+oo,t) -0(-oo,t). (9) 
-oo J — oo 

We see that topological properties are related to large 9 fluctuations and, therefore, the 
constant p approximation should not interfere in such properties. 



In order to explicitly confirm the fact that the soliton operators introduced in 



,13, 



are 



indeed creation operators of quantum solitons in the phase of the field 0, let us explicitly 
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evaluate the commutation relation between the quantum soliton creation operator and the 
topological charge. 

Observe that this is nothing but the Mandelstam creation operator of quantum solitons 
in the SG model [9(, as it should. Then, using canonical commutation relations we readily 
find 

f 27T /*°° f°° "1 

= j -i— / dx'd x ,J d£ x [6{x\t) , 7r e (6, t)]W 

Eq. (jTUJ) implies that the operator /x creates eigenstates of the topological charge Q, with 
eigenvalue 2n/N, thus proving that the quantum solitons occurring in the theory described 
by ([I]) are indeed phase solitons. Corre 



have been calculated elsewhere [1 
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ation functions of these quantum soliton excitations 



121 ] . In Section 4, we will show that the relevant 



quantum correlators for the calculation of the conductivity will be the soliton current-current 
correlators. 



III. THE CASE N=2: A MODEL FOR POLYACETYLENE 

In this section, we are going to propose a phenomenological theory for polyacetylene 
that will enable us to compute quantum soliton correlation functions and in particular the 
correlation functions of quantum solitonic current operators. As we shall see the standard 
field theory model for this polymer unfortunately does not allow the application of the 
method of soliton quantization described in the previous section. For this reason, we will 
propose a phenomenological alternative. 



whose field 



Polyacetylene is described by the Su-Schrieffer-Heeger (SSH) model 
theory version is the Takayama-Lin-Liu-Maki (TLM) model [13J , described by the Hamilto- 
nian 

Htlm = J dxV\{x) \-ihv F a 3 d x + A(x)<r 1 ] *,(ar) 

+ (2nhv F \)- 1 J dx A 2 (x)/Q 2 + A 2 (x) . (11) 

In the above expression, ^ s (x) is a two-component Dirac fermion field, associated to the 
7r-electrons and A fx) is a real scalar field - the dimerization field - associated to the lattice 
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degrees of freedom, namely, the phonons. Furthermore, Oi are the Pauli matrices, Vp is the 
Fermi velocity, A is the dimensionless electron-phonon coupling constant, and is the bare 
optical-phonon frequency. 

Integrating over the fermion field in the previous expression, we obtain an effective theory 
for the phonon field A(x), whose potential is a Z(2) symmetric double- well ig. The 

two degenerate minima of this potential correspond to the two degenerate dimerizations of 
irans-polyacetylene. The equivalent calculation has also been performed within the SSH 
model [3j, |4|, also leading to a degenerate double- well effective potential for the dimerization 
variable. 

The double-well potential for the effective dimerization field implies the existence of 
soliton excitations in the theory. It has been shown that these, indeed, are introduced by 
doping polyacetylene with halogen or alkali atoms 

As we have argued in the Introduction, in the regime of intermediate doping, a full 
quantum treatment of the soliton excitations is unavoidable. Hence, one would be naturally 
inclined to use the method of soliton quantization for theories with a Z(N) symmetry, 
described in the previous section, for the case N = 2. Nevertheless, one immediately realizes 
that the method is not applicable for a real field such as the dimerization field A (2). Indeed, 
for a real field the exponent of the soliton operator ([2]) vanishes, making /x trivial and the 
soliton correlator ([3]) no longer makes sense, since we cannot couple the external field 
to a real field. Furthermore, for a real field, 8 = 0, hence the topological current and the 
respective topological charge cannot be defined as in (jSJ) and ([9]). The above method of 
soliton quantization, in the case of a multiplicative symmetry 6[], only applies to complex 
fields. 

In order to conciliate this fact with the knowledge that the effective theory for the real 
A(x)-field is a degenerate double- well with a Z(2) symmetry, we propose a complex extension 
4>(x) of the A(x)-field, governed by the Lagrangian (CQ) with N = M = 2. The corresponding 
potential is 

V(cj>*, 0) = -7 (0* 2 + 2 ) + V (0*0) 2 , (12) 

or, in terms of polar fields, 

V(p, 9) = -2 7 p 2 cos 26 + rjp 4 . (13) 
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As we can see, there are two degenerate minima at 



(PoA)H V (14) 



where the potential has the value V(p ,9 ) = —■y 2 /rj. Applying, then, the constant p 
approximation 

p(x,t)~p = y| (15) 
and adding jr\ we get the following SG potential for the phase field 9, 

K(0) = -^(l-cos20). (16) 

The associated classical solitonic excitations will be 

6{x) = ±2 arctanexp [y/j (x — xq)] , (17) 

where the plus and minus signs correspond respectively to a soliton (9 s (x — xq)) and an 
anti-soliton in the phase of the 0-field. We see that 

lim 9 s (x — x ) = and lim 9 s (x — xq) = n, (18) 



namely, the phase soliton connects the minima of the potential ( |T3l) . when p = \pyjr}. 

The actual potential for the dimerization field A(x) and the Z(2)-symmetric complex 
potential for the field 0, given by (jT3l) . both possess the same topology, related to the 
Z(2)-symmetry. It is therefore reasonable to expect the same topological properties in both 
theories, especially those concerning solitons. We may adjust the parameters in such a way 
that the minima of the real field potential coincide with those of the complex one. 

What we are doing is quite similar to what is done when we use complex functions in 
order to describe the EM field. The physical E and B fields will correspond to the real part 
thereof. The dimerization field A of the TLM model is the real part of our <fi. As we know, 
for the polyacetilene soliton, we have the A-field varying from — A to A in between the 
two minima. 

A SG soliton in theta would have the phase of the <fi field varying from 7r to 0, implying 
(j) would vary between — p and p . Thus, identifying p with A , we can figure out the 
relation between the sine-Gordon solitons of our model and the polyacetylene solitons: the 
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real part of the complex field for the configuration having a SG soliton in its phase will 
be in the same topological class as the 4 -like soliton of polyacetylene. Since for each SG 
soliton there is a soliton in polyacetylene we may identify the SG-soliton current with the 
polymer soliton current. We are going to use, therefore, this complex extension in order to 
study the quantum properties of the soliton excitations by means of the method of soliton 
quantization described in the previous section. We will use, in particular, the SG soliton 
current for calculating the conductivity. 



IV. THE QUANTUM SOLITON CURRENT CORRELATOR AND CONDUCTIV- 
ITY 



In this section, we are going to obtain an exact series expression for the soliton dc- 
conductiviy in our model for polyacetylene. For this purpose, the starting point is the 



well-known Kubo formula 



13, 



of = lim lim -Im [{f J j ) ret (u, k)] , (19) 

where ( J 1 Ji) ret (u, k) is the retarded, Minkowski space, correlation function of the spatial 
component of the soliton current operator J M given by ([8]). 

We want to evaluate the above current- current correlator within our field theory model 
for polyacetylene. The strategy will be to derive a generating functional for such current 
correlators in our field theory model. For this purpose, we introduce the identity 

1 = J VJ» 5 [J 11 - e" v d u 9] 

in the Euclidian vacuum functional associated to the Lagrangian (J7J), obtaining 

Z = Z l J VJ»V\V6 exp |- J- J d 2 z [($8^6 

-2 7P ^ cos N6 - ihv s (J" - e^d v 6) Aj } , (21) 

where we used expression (|SJ) for the soliton current. Notice that in the expression above, 
we no longer make h = c = 1. Actually, since the 0-field theory replaces the effective theory 
for the phonon field A in the TLM model, we substitute c for v$, the speed of sound in the 
polymer. 
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We now integrate over 9 and A M , thereby obtaining the partition function expressed as 
the functional integral of the exponential of an effective J M action [3]. The 9 integral may 
be done by the usual expansion in powers of the cosine term [ljj] , or equivalently, in powers 
of Pq . The resulting functional integrals, both in 9 and A M are quadratic and the final result 



is 



Z = J\f 



jv \ 2m 

oo I lEa_ 

\ fit's 

m=0 



[mi 



2m 

}\d 2 : 

i=i 

x J*{2 

2m 



VJ" vx]){ -- I J' 2 :,! 2 :' 



+ / d 2 z 



hv s 

iN \^ a d a G(zi - z 



i=l 




(22) 



where Aj = 1 for 1 < i < m and Aj = —1 for n + 1 < i < 2m and G(z) is the Euclidian Green 
function of the two-dimensional (2D) free massless scalar theory, which appears naturally 
since the expansion in powers of the cosine term is an expansion around such theory. 

We note at this point that, should we integrate the above expression over J M , we would 
obtain the usual Coulomb gas representation for the vacuum functional of the SG theory 
[20! ] . Conversely, the expression for the generating functional of current correlators can be 
obtained by the usual procedure of adding a linear coupling with a source in the exponent 
of the integrand in the previous espression, namely, 



Z[K,]=AfJ2 



w \ 2m 
2£o_ 

hv s 



m=0 



(ml)' 



2m 

\\d 2 , 

1=1 

xJ»{z) 



1 



VJ 11 exp<! -- / r/ 2 :r/- : 



21 



2p| 

hv<. 



8\z-z') 



J u (z> 



+ / d 2 z 



2rn 



iNj2^ a 9 a G(zi-z) + 



i=l 



(23) 



Z [K^] is the desired generating functional of J M correlators. Indeed, we have 



(hv s ) 2 S 2 Z[K, 



Z SK^x) 6K U 



(24) 



K u =0 
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Now, integrating ([23]) in J M , we get 



oo I J£o_ 

hv s 



2 in 



m=0 



2m 

IP 



i=l 



x exp 



d 2 zd 2 z' 



2m 



«iV A^ Q d a G(i; - z) + 



i=l 



2p§ 



iN ^^%'G(z j -z') + ^ s 



2m 



3=1 



oo (ifL 



2 m 



x exp 



x exp 
1 



N 2 hv s 



m=0 
2m 2m 

i=l 3=1 



2m 

if 



2m 



iN 



d 2 zd 2 z' K^z)5\z- z')K^z') 



Zp o i=l 



Evaluating the functional derivatives in (1241 . we obtain 



(J"(f)^(y)> = Z- 1 J] 



2m 



m! 



x exp 
kv 



N 2 hv s 



S 5^5 2 {x-y) 



2PI 



m=0 
2m 2m 

i=l 3=1 

iV 2 (^ 5 ) 2 2 "' 



2m 
i=l 



2^ 



4p^ 



2 m 



x 



£ Aj (^(f)^) - G(* + (x - y)) 



i=l 



where we have made the shift of variable z^j) — > Zjy) — x. 
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We now perform the Fourier transform in the variable x = ^ — y, arriving at 

oo ( 1£sl \ „ 2m 



m=0 V '' J i=l 



j2 f- 2m 2m 



x exp \ S Ai 5Z A ^ G( ^ ~ % 



v , Tws xl a, N 2 (hv s ) 2 , . 

2,;i V^ 2 - /,•"/,• 



X 



EV-^^^j)), (27) 
where = (fc, oS) and we used the fact that 

(Siiwt.2 _ 
p J • (28) 

From ( l2Tj) we can get ( J 1 Ji) ret (u), k), by following the prescription given in [17| . for the 
retarded Green function, which includes the change of variables (recalling that z = (z, r) ) 

ir — > —vst, iu) —> h iS, (29) 

v s 

and the limit S — > 0. The soliton conductivity, then, is given by 

<j s = lim lim lim — Im [(JJ) ret (u, k)) . (30) 

5^0 lu^O k^O UJ 

Taking the above limits, we obtain, after some algebra 

r2/fc \2 v I 7P ° I - 2m 



P0 m=l V J i=l 



x exp ^ — zJ^ S Ai S a j' g (^ ~ ^ 



4p^ 

1 v i=1 J=1 



2m ( 2m 



X 



. i=i J I t=i 



The non-transverse part of (|27p is a non-physical "zero-point" term, which must be 
subtracted from the current correlator. Anyway it would not contribute to the conductivity 
because it is real. 
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The temperature dependence of the soliton conductivity may now be obtained by the 
usual methodology, through which we are led to the version of (l3Tj) having finite integration 
regions < Tj < (5 (/3 = llvs/hsT) in the Euclidian time, 



af(T)=ia s (T) 



-iZ- l N 2 (hv s ) 
x exp 



m=l 



/ at \ 2m 
( 1£<L \ r h "S 

\hv s 



N 2 hv s 



ml) 2 j 

2m 2m 



2in 



\\ dn dzi 



i=l 



2m 



i=\ j=l 
2m 



(32) 



i=i 



In the above expression, the free massless Green function has been replaced by the corre- 
sponding function at a finite temperature T, namely, Gt{z) {z = (z,r)). This is a natural 
consequence of the frequency quantization in the presence of a finite interval for r. 

The thermal Euclidian Green function of the 2D free massless scalar theory in coordinate 
space has been evaluated in 2J| and is given by 



G T (z) = --In 



VIP 2 



it 



, , 2nk B T 
cosh — z 



'2nk B T 

COS I — : T 



, -i \ , ■ i . • ( 33 ) 

hvs J V hvs 

Furthermore, the one-dimensional electrical conductivity is related to the quantum soliton 
conductivity presented in ( |32l) by 



<7 P 



e 2 v s 

h 



(34) 



In order to obtain a result that could be compared with experimental data, i.e., the three- 
dimensional electrical conductivity, we must divide the above expression by the cross-section 



? 2 

area of the polyacetylene fibers, namely, A ~ 7T.10 4 A 22]. Then, we have 



a 



A' 



(35) 



In what follows, we will explicitly demonstrate that Eq. (132|) yields an exactly vanishing 
quantum soliton conductivity. First of all, let us change our notation by defining 



K + = ( z t> r i + )> for l <i<m; 



for m + 1 < i < 2m. 



(36) 
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Using the above notation and making the change of variables 
may then rewrite Eq. fl32l as 



af(T)=ia s (T) 



-iZ- 1 N 2 (hv s ) 



x i mi 



where 



E 



2m+l 



V 5 4p4 ^ ( m !)2 ^TrfcflT 

u m=l v ' x 

n y rfr ^ / ^ J s ™( r i + ' • • • ' r i~' • • • ' r * 



(37) 



1=1 m (^l j • • • ) r m i r l ) 



ew-t) n 



i=l 

r2j 



«> — oo J — oo 



¥ 2 o 



i,3=l 



-G T (z+ - zf) - G T (zr - ?+)]} 

m 

.i=i 



in which, after re-scaling r, 



G ^) = -- ln 1 — — 



47T 



7T^ 



, , 2nk B T 
cosh ( — z | — COST 

hv s 



(38) 



(39) 



We will now show that the quantum soliton conductivity vanishes exactly For this, we 
make the change of variables ^ — > 2n — in the r-integrals in fl37|) . Since 



t (27r - r x + , . . . , 27T - r+; 2vr - r x , . . . , 2tt - r m ) 

= — "m( r l j • • • j T'mi T"i > • • • > r m )) 



(40) 



the announced result immediately follows. 



V. CONCLUSION 



The application of a general method of soliton quantization, based on order-disorder 
duality, to a Z(2) symmetric complex extension of the effective field theory for the dimer- 
ization field of the TLM model for polyacetylene has yielded an exactly vanishing result for 
the quantum soliton conductivity This strongly suggests that dynamic solitons are not the 
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charge carriers in polyacetylene in the intermediate doping regime. The natural candidates 
are polarons. However, as it happens in the case of solitons, which were studied in the 
present work, a full quantum treatment is required in order to derive a reliable expression 
for the polaron conductivity as a function of the temperature. We are presently investigating 
the quantum polaronic conductivity in this system. 

This work has been supported in part by CNPq and FAPERJ. LM and AAS were sup- 
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